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Abstract 

We present extended versions and give detailed proofs of results concerning percolation 
(using various sets of two-replica bond occupation variables) in Sherrington-Kirkpatrick spin 
glasses (with zero external field) that were first given in an earlier paper by the same authors. 
We also explain how ultrametricity is manifested by the densities of large percolating clusters. 
Our main theorems concern the connection between these densities and the usual spin overlap 
distribution. Their corollaries are that the ordered spin glass phase is characterized by a unique 
percolating cluster of maximal density (normally coexisting with a second cluster of nonzero 
but lower density). The proofs involve comparison inequalities between SK multireplica bond 
occupation variables and the independent variables of standard Erdos-Renyi random graphs. 
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1 Introduction 



In Ising ferromagnets (with no external field), it is well known that the ordered (broken symmetry) 
phase manifests itself within the associated Fortuin-Kasteleyn (FK) random cluster representa- 
tion yj by the occurrence of a single positive density percolating cluster (see 121). In a recent 
paper (3], we investigated the nature of spin glass ordering within the FK and other graphical rep- 
resentations and concluded that the percolation signature of the spin glass phase is the presence of 
a single two-replica percolating network of maximal density, which typically coexists with a sec- 
ond percolating network of lower density. The evidence presented in that paper for this conclusion 
was two-fold: suggestive numerical results in the case of the three-dimensional Edwards-Anderson 
(EA) spin glass [|4j| and rigorous results for the Sherrington-Kirkpatrick (SK) spin glass 

In this paper, we expand on those results for the SK model in several ways. First, we give much 
more detailed proofs, both for two-replica FK (TRFK) percolation and for the different percolation 
of "blue" bonds in the two-replica graphical representation studied earlier by Chayes, Machta and 
Redner [[6113 (CMR). Second, we go beyond the ±J SK model (as treated in [3J) to handle quite 
general choices of the underlying distribution p for the individual coupling variables, including the 
usual Gaussian case. Third, we organize the results (see in particular Theorems [T] and [2l) in such 
a way as to separate out (see Theorems |4] and [5]) those properties of the overlap distribution for 
the supercritical SK model that are needed to prove related properties about percolation structure. 
Such a separation is called for because many properties of the overlap distribution that are believed 
to be valid based on the Parisi ansatz (see |^) for the SK model have not yet been rigorously 
proved. 

Another way we expand on the results of the previous paper is to present (in Section [3]) an 
analyis of the percolation signature of ultrametricity in the SK model, which is expected to occur, 
based on the Parisi analysis, but has not yet been proved rigorously. That is, we describe (see Theo- 
rems [6l|7] and (8]) how the percolation cluster structure of mutliple networks of differing densities in 
the context of three replicas would exhibit ultrametricity. We note that as a spinoff of Theorem [H 
we have (at least in the SK model — numerical investigations for the EA model have not yet been 
done) a third graphical percolation signature of the spin glass transition beyond the two analyzed in 
our earlier work — namely one involving uniqueness of the maximum density percolating network 
(one out of four clusters) in a three replica mixed CMR-FK representation. 

In addition to these extensions of our earlier results and as we note in a remark at the end of 
this introductory section, the technical machinery we develop in this paper can be used to obtain 
other results for the SK model, such as an analysis of large cluster densities at or near the critical 
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point. Before getting to that, we first give an outline of the other sections of the paper. 

In Section [2] we describe the SK models and the CMR and TRFK percolation occupation vari- 
ables we will be dealing with throughout. We then present our main results, starting with Theo- 
rems [Hand [21 which relate, in the limit ^ oo, the densities of the largest percolation clusters to 
the overlap distribution for the CMR and TRFK representations respectively. Then, after stating a 
known basic result (Theorem [3]) about the vanishing of the overlap in the subcritical (and critical) 
SK model, we present Theorems [4] and [51 which give respectively the CMR and TRFK percolation 
signatures of the SK phase transition, under various assumptions about the SK overlap distribu- 
tion. Results relating percolation structure and ultrametricity in the SK model are presented in 
Section [3l Then in Section |4l we present all the proofs. That section begins with three lemmas 
that are the technical heart of the paper and explain why one can compare, via two-sided stochas- 
tic domination inequalties, SK percolation occupation variables to the independent variables of 
Erdos-Renyi random graphs [9|. A key feature of these comparison results is that they are done 
only after conditioning on the values of the spin variables in all the replicas being considered. 
This feature helps explains why the size of the overlap is crucial — because it determines the sizes 
of the various "sectors" of vertices (e.g., those where the spins in two replicas agree or those where 
they disagree) within which the comparisons can be made. 

Remark 1.1 The first part of Theorem\5\says that for [3 <1, the size of the largest TRFK doubly 
occupied cluster is o{N) or equivalently that its density D( = o{l) as N ^ oo. But the proof 
( see the TRFK part of Lemma \4.2\) combined with known results about G{N, p^), the Erdos-Renyi 
random graph with N vertices and independent edge occupation probability pn, implies quite a 
bit more — that dJ^^^ is 0{\ogN/N) for (3 < 1 and 0{N'^'^/N) for [3 = 1. Even more is 
implied in the critical case. E.g., in the critical scaling window, where (3 = (3^ = 1 + X/N^^"^, 
the largest clusters, of size proportional to N'^^^, behave exactly like those occuring in a pair of 
independent random graphs ( see, e.g., / [70] [771 [721 [731/ ) — i.e., as N ^ oo, the limiting distribution 

is the same as that obtained by taking two independent copies 
of G{N/2,2{Pn)^/N), combining the sizes of the largest clusters in the two copies, then rank 
ordering them and dividing by N'^/'^. One can also show that for [3 > 1, the size of the third largest 
TRFK cluster behaves like that of the second largest cluster in a single copy of the supercritical 
Erdos-Renyi random graph — i.e., 0(\ogN) /[9l/. But that derivation requires a strengthened 
version of Lemma 14. 21 and further arguments, which will not be presented in this paper. 
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2 Main Results 



Before stating the main results, we specify the random variables we will be dealing with. For 
specificity, we choose a specific probabilistic coupling so that even though we deal with two dif- 
ferent graphical representations, and a range of inverse temperatures (3, we define all our random 
variables for the system of positive integer size on a single probability space. The corresponding 
probability measure will be denoted Pjv (with P denoting probability more generically). 

For each A^, we have three types of random variables: real- valued couplings { Jjj}i<i<j<Ar, 
Ising ±l-valued spins {(Ji}i<i<N and {Ti}i<i<i\i for each of two replicas, and a variety of percola- 
tion {0, l}-valued bond occupation variables which we will define below. These random variables 
and their joint distributions depend on both N and /? (although we have suppressed that in our no- 
tation), but to define them, we rely on other sets of real- valued random variables not depending on 
N or (3: {-ft'ij}i<i<j<oo and {C//j}i<j<j<oo for each replica indexed by £ = 1 or 2. (In later sections, 
we will consider more than two replicas.) Each of these sets is an i.i.d. family and the different sets 
are mutually independent. The f//^'s are independent mean one exponentials and will be used to 
define the bond occupation variables (conditionally on the couplings and spins). The i^'y 's, which 
determine the 's for given N (and (3) by Jij = Kij/ \/N, have as their common distribution a 
probability measure p on the real line about which we make the following assumptions: p is even 
(dp{x) — dp{—x)) with no atom at the origin (p({0}) = 0), variance one {j^^ x^dp{x) — 1) and 
a finite moment generating function {j^^ e^^dp{x) < oo for all real t). The two most common 
choices are the Gaussian (where p is a mean zero, variance one normal distribution) and the ± J 
(where p is {6i + S-i)/2) spin glasses. 

For a given N and /3, we have already defined the couplings Jy. The conditional distribution, 
given the couplings, of the spin variables cr, r for the two replicas, is that of an independent sample 
from the Gibbs distribution; i.e.. 



const X exp 



l<i<j<N 



(1) 



It remains to define the percolation bond occupation variables of interest, given the couplings 
and the spins. Two of these are the FK (random cluster) variables — one set for each replica; we 
will denote these nf^ for i = 1 (corresponding to the first (a) replica) and i = 2 (corresponding to 
the second (r) replica). These may be constructed as follows. For a given i < j, if the bond {i, j} 
is unsatisfied in the first replica — i.e., if JijCTiCTj < 0, then set njj = 0; if the bond is satisfied, 
then set n,]^ = 1 if Ulj < 2P\Jij\ (i.e., with probability 1 — cxp (— 2/3| Jjj|)) and otherwise set it 
to zero. Define n? similarly using the second (r) replica. We will be particularly interested in the 
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percolation properties of the variables rijj = T^jj^fj that describe doubly FK-occupied bonds. We 
will use the acronym TRFK (for Two Replica FK) to denote various quantities built out of these 
variables. 

There is another two-replica graphical representation, introduced by Chayes, Machta and Red- 
ner EC! (which we will denote by CMR) that we will also consider. This representation in general 
is described in terms of three types of bonds which may be thought of as those that are colored 
blue or red or else are uncolored. One way of defining the blue bonds, whose occupation variables 
we will denote by bjj, is that bjj = 1 if {«, j} is satisfied in both replicas and also either nj^ = 1 
or = 1 or both (which occurs with probability 1 — exp ( — 4/? | Jj^ | ) ) ; otherwise bj^ = . We will 
be interested in the percolation properties of the blue bonds. Although we will not be using them 
in this paper, we note that a bond {i, j} is colored red if and only if aiajTiTj = —1 (or equivalently 
{i,j} is satisfied in exactly one of the two replicas) and also that satisfied bond is FK-occupied 
(which occurs with probability 1 — exp (— 2/9| Jjj|)). 

A key role in the theory of spin glasses, and this will also be the case for their percolation 
properties, is played by the Parisi (spin) overlap. For a given N and (3, this overlap is the random 
variable, 

Q = Q{N,(3) = J2 (2) 

l<i<N 

Closely related to the overlap are the densities (i.e., the fractions of sites out of A^) Da = Da{N, (3) 
and Dd = Dd{N,P) of the collections of sites where the spins of the two replicas respectively 
agree and disagree with each other. Since Q = Da — Dd and Da + Dd = 1, one can express 

Dmax = max{Da, Dd} and Dmin = m.m{Da, Dd} as Dmax = [1 + |Q|]/2 and Dmin = [1 - IQIl/2. 

It should be clear from our definitions of the various bond occcupation variables that if one of 
i, j is in the collection of agree sites and the other is in the collection of disagree sites, then the bond 
{i,j} is satisfied in exactly one of the two replicas and so {i,j} can neither be a TRFK occupied 
bond nor a CMR blue bond. So percolation (i.e., occurrence of giant clusters containing order A^ 
of the sites) can only occur separately within the agree or within the disagree collections of sites. 
Our results concern when this happens and its connection with the spin glass phase transition in 
the SK model via the overlap random variable Q. 

We will first state two general theorems relating the occurence of giant clusters to the behavior 
of Q and then state a number of corollaries. The corollaries depend for their applicability on results 
about the nature of Q in the SK model, some of which have been and some of which have not yet 
been derived rigorously. The first theorem concerns CMR percolation. We denote the density of 
the fc'th largest CMR blue cluster by D^^^{N, {3). 
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Theorem 1 In the CMR representation, for any < /5 < oo, the following three sequences of 
random variables tend to zero in probability, i.e., the Y* is[ -probability that the absolute value of the 
random variable is greater than e tends to zero as N ^ oo for any e > 0. 

(3) 

DP^(iV,/?)-[i^M^]^0. (4) 

D^M^{N,P)^0. (5) 

To state the next theorem, we define 6{c) for c E [0, oo) to be the order parameter for mean- 
field percolation — i.e., the asymptotic density (fraction of sites out of as oo) of the 
largest cluster in G{N, c/N), the Erdos-Renyi random graph with occupation probability c/N 
independently for each edge in the complete graph of sites [9j. It is a standard fact (9] that 6{c) 
is zero for < c < 1 and for c > 1 is the strictly positive solution of 

6 = 1- e-"^ . (6) 

Theorem 2 In the TRFK representation, for any < /3 < oo, the following limits (in probability) 
are valid as N oo for D^^^^, the density of the k'th largest TRFK doubly occupied cluster 

DP^^(iV.3)-0(2fl^l±MiMa)|i±M^]^0. (7) 

dWHn,0) - ei20^' - l^f -^'l)!' - "^f •'^"l ^ 0. (8) 

D™'''^(N,I}) ^0. (9) 

We now denote by P{N,P) the probability distribution of the overlap Q(N,P). This is the 
Parisi overlap distribution, averaged over the disorder variables /C = {-ft'ij}i<i<j<oo- The unaver- 
aged overlap distribution requires conditioning on /C. So, for example, we have, for g G [—1,1], 

P(Ar,/3)([-l,g]) = Av[P;v(Q(iV,/3) <g|/C)], (10) 

where Av denotes the average over the disorder distribution of JC. 

The quantity E{Q{N,py\}C) is closely related to J2i<i<j<N[Eiaiaj \ }C)f (see, e.g., il4J . 
Lemma 2.2.10) which in turn is closely related to the derivative of the finite volume free energy 
(see, e.g., IBl, Prop. 4.1). It then follows that E{Q{N, p^) ^ as N ^ oo first for /3 < 1 ([15], 
Prop. 2.1) and then (using results of [[H [HI [H and of [[13 |20l| — see also [21J) also for (3 = 1. 
This implies the following theorem, one of the basic facts in the mathematical theory of the SK 
model. 
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Theorem 3 For (3 <1, Q{N, P) ^ (in probability) as N ^ oo, or equivalently P{N, (5) 6o. 

The situation regarding rigorous results about the nonvanishing of Q(N,P) as ^ oo for 
P > lis less clean. For example, using results from the references cited just before Theorem [3l 
it follows that E{Q{N, P)"^) has a limit as ^ oo for all /5, related by a simple identity to 
the derivative at that P of the infinite-volume free energy, V{P), given by the Parisi variational 
formula, V{P) = inf^ V{m, P), where the inf is over distribution functions m{q) with q E [0, 1] 
— see [[T9l . Furthermore, since for P > 1, V{P) is strictly below the "annealed" free energy II22I 
(which equals V{5q, P), where 5q is the distribution function for the unit point mass at g = 0), it 
follows (see [21]) from Lipschitz continuity of P(m, P) in m |fT7ll20| that 



However, it seems that it is not yet proved in general that Q{N, P) has a unique limit (in distri- 
bution) nor very much about the precise nature of any limit. In order to explain the corollaries 
of our main theorems without getting bogged down in these unresolved questions about the SK 
model, we will list various properties which are expected to be valid for (at least some values of) 
P > 1 and then use those as assumptions in our corollaries. Some related comments are given in 
Remark 12.21 below. 

Possible Behaviors of the Supercritical Overlap. For P > 1, P{N, P) converges as ^ oo 
to some Pj3 with the following properties. 

• Property PI: P/3({0}) = 0. 

• Property P2: P/3({-l, +1}) = 0. 

. Property P3: P^([-l, -1 + (l//?^)] u [1 - (l//?^), 1]) = 0. 

If one defines QEAiP), the Edwards-Anderson order parameter (for the SK model), to be the supre- 
mum of the support of Pp, and one assumes that has point masses at ±qEAiP), then Properties 
P2 and P3 reduce respectively to QEAiP) < 1 and QEAiP) < 1 — {l/P'^). Weaker versions of the 
three properties that do not require existence of a limit for P(A^, P) as N ^ oo are as follows. 

• Property PI': 



lim E{Q{N,pf) > forall /? > 1. 



(11) 



limlimsupP7v(|(5(iV,/3)| < e) = 0. 



(12) 



• Property P2': 



limlimsupP7v(|Q(A^,/?)| > 1 - e) = 0. 



ej.0 N^oo 



(13) 
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• Property P3' 

limlimsupP7v(|Q(Ar,/5)| > I - {I / fS"^) - e) = 0. (14) 

We will state the next two theorems in a somewhat informal manner and then provide a more 
precise meaning in Remark IXTI below. 

Theorem 4 (Corollary to Theorem [1]) In the CMR representation, for any < /3 < 1, there 
are exactly two giant blue clusters, each of (asymptotic) density 1/2. For 1 < /3 < oo, there are 
either one or two giant blue clusters, whose densities add to 1; there is a unique one of (maximum) 
density in (1/2, 1] providing Property PI (or PI') is valid and there is another one of smaller 
density in (0, 1/2) providing Property P2 (or P2') is valid. 

Theorem 5 (Corollary to Theorem lU In the TRFK representation, there are no giant doubly 
occupied clusters for j3 < 1. For 1 < /3 < oo, there are either one or two giant doubly occupied 
clusters with a unique one of maximum density providing Property PI (or PV) is valid and another 
one of smaller (but nonzero) density providing Property P3 (or P3') is valid. 

Remark 2.1 For0<f3<l, Theorem^states that {D^^^{N, (3), D^^^{N, (3), Df^^{N, f3)) 
converges (in probability or equivalently in distribution) to (1/2, 1/2,0) while Theorem\5\ states 
that the corresponding triple of largest TRFK cluster densities converges to (0, 0, 0). A precise 
statement of the results for [3 G (1, C)o) is a bit messier because it has not been proved that there is 
a single limit in distribution of these cluster densities, although since the densities are all bounded 
(in [0, 1]) random variables, there is compactness with limits along subsequences of N's. For 
example, in the CMR case, assuming Properties PI' and P2', the precise statement is that any limit 
in distribution of the triplet of densities is supported on {(1/2 + a, 1/2 — a, 0) : a G (0, 1/2)}. 
Precise statements for the other cases treated in the two theorems are analogous. 

Remark 2.2 Although Property PI' does not seem to have yet been rigorously proved (for any 
P > I), a weaker property does follow from ([77]). Namely, that for all (3 > 1, the limit in ([72]) 
is strictly less than one. Weakened versions of portions of Theorems^and\5\for (3 > 1 follow — 
e.g., any limit in distribution of the triplet of densities in Remark \2J\ must assign strictly positive 
probability to {(1/2 + a, 1/2 - a, 0) : a G (0, 1/2]}. 

3 Ultrametricity and Percolation 

In this section, in order to discuss ultrametricity, which is expected to occur in the supercritical 
SK model (see we consider three replicas, whose spin variables are denoted {erf} for i = 
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1,2,3. We denote by nfj the FK occupation variables for replica £ and by bf™ the CMR blue bond 



occupation variables for the pair of replicas £, m. Thus b^^ corresponds in our previous notation 



to hij. We also denote by Q^*" = Q^^{N, (3) the overlap defined in Q, but with a, r replaced by 

Let us denote by P^{N,(3) the distribution of the triple of overlaps {Q^^ , Q^^ , Q'^^) ■ Ultra- 
metricity concerns the nature of the limits as — > oo of P^(N, as follows, where we define 

^ultra ^ {(^'2/5-2) : \x\ = \y\ < k|or|a;| = \z\ < \y\or\y\ = \z\ < \x\} . (15) 



Possible Ultrametric Behaviors of the Supercritical Overlap. For f3 > 1, P^{N,(3) con- 
ges to some Pg as oo with 

. Property P4:P|(R3^^ J = 1. 



We will generally replace this property by a weakened version, P4', in which it is not assumed that 
there is a single limit P| as ^ oo but rather the same property is assumed for every subsequence 
limit. There is another property that simplifies various of our statements about how ultrametricity 
is manifested in the sizes of various percolation clusters. This property, which, like ultrametricity, 
is expected to be valid in the supercritical SK model (see [[231, where this property is discussed and 
also numerically tested in the three-dimensional EA model) is the following. 

• Property P5: Pp{{{x,y, z) : xyz > 0}) = 1. 

Again we will use a weaker version P5' in which it is not assumed that there is a single limit P| as 

A^ oo. 

One formulation of ultrametricity using percolation clusters is the next theorem, an immediate 
corollary of TheoremfH in which we denote by Dj™ = Dj'^(N, (3) the density of sites in Cj™, the 
j'th largest cluster formed by the bonds {i, j} with bf™ = 1 (i.e., the j'th largest CMR blue cluster 
for the pair {£, m} of replicas). Note that coincides in our previous notation with D^^^. 

Theorem 6 (Corollary to Theorem [1]) For 1 < (3 < oo, assuming Property P4', any subsequence 
limit in distrbution as N oo of the triple (Z^P — D^"^, — D^^, Df" — D^) is supported on 

ultra' 

In our next two theorems, instead of looking at differences of densities, we express ultrametric- 
ity directly in terms of densities themselves. This is perhaps more interesting because rather than 
having three density differences, there will he four densities. We begin in Theorem |7] with a fully 
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CMR point of view with four natural non-empty intersections of CMR blue clusters. Then Theo- 
rem[8]mixes CMR and FK occupation variables to yield (four) other natural clusters. 

There are a number of ways in which the four sets of sites in our fully CMR perspective can 
be defined, which turn out to be equivalent (for large N). One definition is as follows. For a, a' 
each taken to be either the letter a (for agree) or the letter d (for disagree), define A^^' (A^, (3) to be 
the set of sites i E {1, . . . , A^} where aj agrees (for a = a) or disagrees (for a = d) with af and 
aj agrees (for a' = a) or disagrees (for a' = d) with af; also denote by Daa'iN, 13) the density 
of sites (i.e., the fraction of A^) in Aaa'iN,P). Then denote by C^^, the largest cluster (thought 
of as the collection of its sites) formed within Aaa' by the b^"^ = 1 blue bonds. Finally, define 
Caa' = Cll, n Cll,, D^^^ to be the density of sites (fraction of A^) in Ca,a' and D^^^{N, (3) to 
be the vector of four densities {D™^, ^S^^, ^£^^, ^£^^)- To state the next theorem, let 
R'^'^ denote {(xi, X2, Xs, X4) : eachxj > 0} and define 

^ultra ^ {(3^1, 2^2, 3:3, 0:4) G R^'+ : X(i) > X(2) > a^(3) = a;(4)}, (16) 
where a;(2), X(3), X(4) are the rank ordered values of xi, 0:2, xs, X4. 

Theorem? For < p < 1, D^^^{N,P) (1/4,1/4,1/4,1/4) (in probability) as N ^ 
00. For 1 < /5 < 00 and assuming Properties PI', P2', P4',P5', any limit in distribution of 
jjCMR(^j^^ /3) is supported on R^'^+^^ . 

Remark 3.1 The equilateral triangle case where \Q^'^\ = \Q^^\ = \Q^^\ = Qu corresponds to 

= (1 + 3g„)/4 and X(2) = X(3) = X(4) = (1 — qu)/^- The alternative isosceles triangle 
case with, say, \Q^'^\ = Qu > IQ^'^l = IQ'^'^l = Qe corresponds to a;(i) = (1 + g„ + 2g^)/4, 
X(2) = (l + Qu- 2qe)/Aand x^s) = X(4) = (1 - g„)/4. 

In the next theorem, we consider C*^/ defined as the largest cluster within Aq,q,/(A^, P) formed 
by bonds {i,j} with b||nfj = 1 — i.e., bonds that are simultaneously CMR blue for the first two 
replicas and FK-occupied for the third replica. Let D*(N, (3) denote the corresponding vector of 
four densities. As in the previous theorem, we note that there are alternative, but equivalent for 
large A^, definitions of the clusters and densities (e.g., as the four largest clusters formed by bonds 
{i, j} with bj|nf^- = 1 in all of {1, ... , A^} without a priori restriction to Kaa' (N, 13)). 

Theorem 8 For < /3 < 1, Z)*(A^, {3) (0, 0, 0, 0) (in probability) as N ^ 00. For 1 < (3 < 00 
and assuming Properties PI', P3', P4', P5', any limit in distribution of D*{N, [3) is supported on 

ultra 
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Remark 3.2 Here the limiting densities are of the form = 0{A[3'^di^j)) d(^j) where d(i) = (1 + 
qu + 2g^)/4, d(2) = (1 + — 2g^)/4 and (i(3) = (i(4) = (1 — Q'm)/4, with = qufor the equilateral 
triangle case. 

4 Proofs 

Before giving the proofs of our main results, we present several key lemmas which are the technical 
heart of our proofs. We will use the notation >> and << to denote stochastic domination (in 
the FKG sense) for either families of random variables or their distributions. E.g., for the m- 
tuples X = (Xi, . . . , Xm) and Y = (Fi, . . . , Ym) we write X « Y or Y » X to mean 
that 'E{h{Xi, . . . , Xm)) < E(/i(Fi, . . . , Y^)) for every coordinatewise increasing function h (for 
which the two expectations exist). 

Our key lemmas concern stochastic domination inequalities in the fc-replica setting comparing 
conditional distributions of the couplings {Kij} or related bond occupation variables, when the 
spins a^, . . . ,a'' are fixed, to product measures. These allow us to approximate percolation vari- 
ables in the SK model by the independent variables of Erdos-Renyi random graphs when ^ oo. 
Given probability measures on R for 1 < i < j < N, we denote by Prod7v({'^j/; }) the cor- 
responding product measure on R^(^~i)/2. We also will denote by p[y] the probability measure 
defined by dp[-f]{x) = e^^rfp(x)//+^ e^^'rfp(a;'). 

Lemma 4.1 Fix N, P, k and let /i^^ denote the conditional distribution, given the spins cr^ = 
{a} : I < i < N), . . . , c^^ of {K[- = eijKij}i^i^j^N (where the Si/s are any given ±1 values). 
Then 

Prod^{{p[^',f]}) « p%^^ « Prodr,{{pb',-^]}) , dV) 

where 

lit = + • • • + ^M) ± ^] • (18) 

Proof. Define the partition function 

Z^,p = Z^,^{{Ki^})=J2^M4^ E (^cr.cr,)^;.) , (19) 

where J2a denotes the sum over all 2^ choices of cTj = ±1 for 1 < i < A^. Thus the normalization 
constant in Equation ([U for two replicas is {Z^,/})^'^ and the A; -replica marginal distribution for 
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{Klj} is as follows, where ProdAr({p}) denotes ProdAr({z/y}) with Uij = p, (p^j = eij{a}aj + . . . + 
(T^aj) and Cn,i3 is a normalization constant depending on the ^f^'s but not on the K'^j's. 

dfl%^p = C^AZNAiK,}))-' exp{^ Y: 0{,.ir;.)Prod;v({p}). (20) 

It is a standard fact about stochastic domination that if /i is a product probability measure 
on and g is an increasing — i.e., coordinatewise nondecreasing — (respectively, decreasing) 
function on R™ (with / gdjl < oo), then the probability measure jig defined as djj,g{xi, . . . , Xm) = 
g{xi, . . . ,Xm)djj,/ J gdjl satisfies jig » jl (respectively, jig « jj). This follows from the fact 
that product measures satisfy the FKG inequalities — i.e., for / and g increasing / fgdjl > 
I fdji J gdjl. 

On the other hand, since each eija^aj = ±1, it follows from (fT9l) that 

with each tp^ = or ±2 and hence each -ip^j > (respectively, each < 0). Thus, as a 
function of the K'^j's, is increasing and (Z^^^)'^ is decreasing while is decreasing and 
(Z^^)~^ is increasing. Combining this with the previous discussion about stochastic domination 
and product measures, we see that 

C^,^e-(^/^)S(-'=^>'^+<^^^-^»'^)prod^({p}) « ji%^^« C+^e"(^/^)S(+'=^»'^+^^^-^»'^)prod^({p}), 

(22) 

which is just Equation (flTI) since the are normalization constants. This completes the proof 
of Lemma 14.1 [ 

The next two lemmas give stochastic domination inequalities for three different sets of occu- 
pation variables in terms of independent Bernoulli percolation variables. The first lemma covers 
the cases of CMR and TRFK variables involving two replicas and the second lemma deals with 
mixed CMR-FK variables involving three replicas. The parameters of the independent Bernoulli 
variables used for upper and lower bounds are denoted p*j^p, where * denotes CMR or TRFK or 3 
(for the mixed CMR-FK case) and [\ denotes u (for upper bound) or £ (for lower bound) and are as 
follows. 

roc 

P*N,p = 9*NAx)dp*^p{x) , (23) 
P*N,p = 9*N,p{x)dp{x) , (24) 

where 
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gJf'H^) = (1 - e~'W^>y , (26) 

g3^^^(x) = (1 - e-^^'^/^)^) (1 - e-2(/3/v^)-), (27) 

and 

pS^^ = Plf^ = P[4/?/v^] , (28) 

P%,p = p[6/9/v^] . (29) 

Lemma 4.2 FixN, j3, cr^, and consider the conditional distributions fi^^^ of{hij = ^ij}i<i<j<N 
andfilf^p o/{nij = n-^nf j}i<j<j<Ar. T/zen 

ProdNip*j6i + (1 - Ki)^o) << P^,/3 << ProdNip*j6i + (1 - P*;/)5o) , (30) 

where p*f = p*j^^p ifi, j are either both in or both in andp*'^ = otherwise. The asymptotic 
behavior as N ^ oo of the parameters appearing in these inequalities is 

PnT'^ = 4|( P WpixW + 0{^)) , (31) 



TRFK,^ 2/32 1 
P^,, •" = ^(l + 0(^)). (32) 

Proof. We begin with some considerations for the general case of k replicas and arbitrary eij 
before specializing to what is used in this lemma. The case k = 3 will be used for the next lemma. 
Let {K.,IJYn = {Kij,Ulj, . . . ,Ufj)N denote random variables (with 1 < i < j < N) whose 
joint distribution is the conditional distribution of {K[j, Ulj, . . . , Ufj)^ given , . . . ,a^. Note that 
the f//j 's are independent of the -ft'j/s and (like the ?7f/s) are independent mean one exponential 
random variables. Let {K^'^^li^ij}], t/i"^)^ = {K\f%ijl . . . , Ulf^'^)^ denote mutually 

independent random variables with _ft'|j^^[7ij] distributed by p[^ij\ and each Ulf^'^ a mean one 
exponential random variable. 
Lemma |4m says that 

{K'''\l'n\)N « iK)% « (^^"^[{7jn])^ (33) 
and it immediately follows that 

{K^^lil'f}], -U'''^)'n « (K, ~U)% « (iri"^[{7S'+}], -U'''^)% ■ (34) 

We now take k = 2, choose Sij = aja^ and note that rijj = nj^n?^ = = bjj unless i,j 
are either both in or both in A^, in which case we have ajaj = crfcTj, '-ffj^ = 4/3/ y/N and 
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"jij = 0. Furthermore nf^ is the indicator of the event that 2(3K-j + {—U-j) > while hij = hjj 
is the indicator of the event that njj + nfj > 1 and K'^j > 0, so that these occupation variables (for 
such are increasing functions of {K', —U)\. Let us now define fiy , bj^ and n|^^'^, b|^^'^ as 
the same increasing functions, respectively, of {K, —U)% and (^('^'^^[{7^^'^}], —U^^^Yj^ providing 
i, j are either both in Aq or both in Ad, and otherwise set these occupation variables to zero. Then, 
as a consequence of (l34l) , we have 

(ni"d,-)2^« (f,)2^« (ni"^'^)^, (35) 



and 



To complete the proof of Lemma l4~2l it remains to obtain the claimed formulas for the various 
Bernoulli occupation parameters. E.g., for the case * = TRFK and t] = m, we have for i,j either 
both in Aa or both in A^^ , 

\l-e-^^^/V^^'')^dp[A(3/y/N]{x), (37) 





as given by Equations (l23l) . (l26l) and (|28|) . We leave the checking of the other three cases (for 
k = 2) and the straightforward derivation of the asymptotic behavior as ^ 00 for all the 
parameters to the reader. This completes the proof of Lemma 1421 

Lemma 4.3 FixN, /3, ct\ a^, and consider the conditional distribution jl]^ ^ o/{bj jnf^ }i<j<j<7v. 
Then equation ^JUi with * = 3 remains valid, where pl'^ = p]^jj if hi (^^^ either both in Aaa or 
both in Aad or both in Ada or both in Add and p^f = otherwise. The asymptotic behavior as 

N 00 ofp]!^p is 

/.', = f + (38) 

Proof. As in the proof of Lemma 14^21 we apply Lemma 1431 in the guise of (|34l ) with eij = ajaj 
again, but this time with k = 3. Now bj^nf^ = ^jfi^ij = unless i, j are either both in Aaa or both 
in Aad or both in Ada or both in Add, in which case we have alaj = afa^ = crfo-j, 7^^-''' = 6/5/ y/N, 
'yfj~ = 0, and ^jf^^^j is an increasing function of {K', — f/)^. The remainder of the proof, which 
closely mimics that of Lemma l4.2[ is straightforward. 
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Proof of Theorem SI We denote by D^^^^ = ^^^^(iV, /3) for a = a or the density (fraction 
of A^) of the j'th largest CMR cluster in A^, and by sf^^ = Sf'^{N,pN) the number of sites in 
the largest cluster of the random graph G{N, p^). Recall that = Da{N, 13) denotes the density 
(fraction of A^) of and that max{L)„, Dd} = [1 + \Q\]/% min{L)„, Dd] = [1 - |<5|]/2. Then 
Lemma l4~2l implies that, conditional on cr\ cr^, 

ND^f\N,(3) » sf^iN D^iN,P),i2Pcp/VN) + Oil/N)) (39) 

with Cp > 0. By separating the cases of small and not so small Da{N,P) and using the above 
stochastic domination combined with the facts that m^^ sf^^ {m, — ^ 0{1) = 1 (in probability) 
when m, mpm oo, and that L'^r^ < D^,it directly follows that - Z^^f^^ and hence 
also that -D^^^ (in probability). This then directly yields © and © and completes the 
proof of Theorem [U 

Proof of Theorem m The proof mimics that of Theorem [T] except that (|39l) is replaced by 

S^^iND^, ^ + OiN-'/^)) » NDjf^^ » S^^iND^, ^ + 0{N-''/^)) , (40) 

where the two terms correcting 2(3'^ /N, while different, are both of order 0(A^~^/^). We then use 

the fact that m^^ Sf^^ {jn,prn) — 9{mpm) when m ^ oo to conclude that dJ^^^{N, [3) — 

9{2p'^Do,{N, P))Do,{N, P)) which yields © and ®. To obtain ®, we need to show that 

Da2 0. But this follows from what we have already showed, from the analogue of (|40l) 

with Z^JP^ replaced by ojf^^ + Z^JP^ and sf^^ replaced by Sf^^ + Sf^^, and by the fact 
thatm i^l^G 

(m, pm) — in probability when m — oo. This completes the proof of Theorem[2l 
Proof of Theoremll Theorem[I]impliesthatany limit in distribution of (Dp^^, D^^^, D^^^) 
coincides with some limit in distribution of (^-^^, 0). The claims of the theorem for < 
P < 1 then follow from Theorem [3] and for 1 < < oo from that theorem and Properties PI' and 
P2'. 

Proof of TheoremlH Theorem[2]implies that any limit in distribution of {dJ^^^, dJ^^, dJ^^) 
coincides with some limit in distribution of {e{f3\l + \Q\)){ ^ ),e{(3\l - \Q\)){ ^ ) , ) . The 
claims of the theorem for < /? < 1 then follow from Theorem [3] and the fact that 9{c) = for 
c < 1 while the claims for 1 < /3 < oo follow from Theorem [3] and Properties PI' and P3'. Note 
that P3' is relevant because 6'(/?2(l - |(5|)) > if and only if /^^(l - |(5|) > 1 or equivalently 
IQI > 1 - (1//32). 

Proof of Theorem HI This theorem is an immediate consequence of Theorem \T\ (which can be 
applied to CMR percolation using any pair of replicas i, m) since — = 
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Proof of Theorem 121 The identity Q^"^ = Da — Da which involves only two out of three replicas 
may be rewritten in terms of three-replica densities as 

Q^^ = Daa + Dad-Dda-DM. (41) 

The corresponding formulas for the other overlaps are 

Q^"" = Daa - Dad + D^a " Dm , (42) 
= Daa - Dad " Dda + Ddd ■ (43) 

Combining these three equations with the identity Daa + Dad + Dda + Ddd = 1, one may solve for 
the Daa''s. to obtain 

Daa = {Q'' + Q'' + Q'' + l)/^, (44) 

Dad={Q'^ -Q'^ -Q^' + 1)/^, (45) 

Dda = i-Q'^ + g^' - g'' + 1)/4 , (46) 

DdA = i-Q'^ - Q'' + + l)/4 . (47) 

In the equilateral triangle case where |g^^| = \Q^^\ = |g^^| and also Q^'^Q^^Q'^^ > (from 
Property P5'), one sees that the corresponding values of the D^^a''^ lie in R-^'i^j.^ with > 
X(2) = = X(4) . In the isosceles triangle case (with Q^^Q^^Q^^ > 0), one is again in R^'+^^, 
but this time with a;(i) > X(2) > a;(3) = a;(4). 

Proof of Theorem m Here we use Lemma 143] to study the densities (fractions of A^) D^^, j (where 
a and a' are a or d) of the j'th largest cluster in A^a' formed by the bonds {ij} with bj^n^^ = 1. 
Similarly to the proof of Theorem [5l we have 

and use this to conclude that D^,^,^^ — 6 {4(3"^ Daa') Daa' and D^^, 2 0. Noting that 
D 6{A[3'^D) is an increasing function of D, the rest of the proof follows as for Theorem|71 
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